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ABSTRACT 

Critical mass flux and axial pressure profile data 
for fluid nitrogen are presented for N - 20, 15, 10, 
and 7 N-sequential-orifice-inlet configurations uni- 
formly spaced at 15.5 cm. These data correlate well 
over a wide range in reduced temperature (0.7 < Tr,o < 
ambient) and reduced pressure (to Pr * 2) and are in 
general agreement with previous studies of one to four 
inlets. ' 

c 

Experimental and theoretical agreement is good for 
liquid and gas critical mass flux, but inconclusive in 
the near-thermodynamic critical regions. 

INTRODUCTION 

The inlets to many flow devices are not the smooth 
entrances often researched. Frequently they contain a 
series of constrictions which may not act independently 
of one another. This class of flows is categorized as 
having sequential inlets. Compressors, pin-finned heat 
exchangers, separators, and labyrinth and step seals 
are examples of inlet configurations consisting of two 
or more such constrictions or, in a stricter sense, 
sequential inlets. The details of the flow dynamics 
and heat transfer in these configurations are in many 
cases not well understood. 

Similarity principles can often be applied to 
assist the designer in cases where the theory is incom- 
plete or no data are available. In Ref. (1) , the prin- 
ciple of similarity as applied to thermophysical pro- 
perties and fluid mechanics (with emphasis on the 
thermodynamic critical region) is shown to qual itativley 
group the experimental results for several physical 
processes, including heat transfer and two-phase choked 
flows. 

In other work (2-5) it was assumed that some form 
of similarity relationship between mass flux and reduced 
pressure was valid. It was then shown analytically and 
experimentally that critical mass flux through well- 
separated sequential inlets for N < 4 could be related 
to flow through a single inlet. The statement appears 
valid even for such diverse inlets as those of the ori- 
fice and Borda types (6). 

These results also show that a similarity exists 
between the flow losses in a rough tube and those in 
sequential inlets. For example, consider a tube of 
length L that is artifically subdivided into N connected 
segments each of length L/N and an N-sequential inlet 
configuration. When the values of the friction parameter 
4fL/D are of order one or more for each tube segment, the 
flow losses through these N connected tubes can be made 
equivalent to the flow losses through the N sequential 
inlets (6,7). 

The primary purpose of this paper is to extend the 
range of available critical (choked) flow and pressure 
profile data of N sequential orifice inlets with N > 4 
for a large range of inlet stagnation conditions. As a 
secondary purpose, existing similarity principles will 
be used and others developed as required to enable de- 
signers to predict the critical mass flow through N se- 
quential inlets and to provide assessments of the pres- 
sure profi les. 


APPARATUS AND INSTRUMENTATION 

The basic components and operations of the blowdown- 
type facility described in Refs. (7) and (8)were modified 
to accommodate the various sequential-orifice-inlet con- 
figurations (Fig. 1). 

A photograph illustrating the N-inlet test configu- 
rations, as installed in the facility, for N - 7, 10, 

15, and 20, is shown in Fig. 2. 

The orifice-type inlets with t/D of 0.5, similar 
to those of Refs. (6) and (7), were designed with spacers 
of 15.24 cm (6.0 inT). These provided fixed spacings of 
32 orifice diameters between apertures, or approximately 
15.5 cm aperture to aperture. A schematic of a segment 
of an N-sequential -orifice-in let configuration illustra- 
ting the pressure tap locations and the inlet geometry 
is presented in Fig. 3; a photograph of the orifice 
proper is shown in Fig. 4. 

These sequential-inlet configurations were fitted 
between inlet and outlet flange adapters to accommodate 
the multiple lengths. Six threaded rods and two plastic 
spacers per segment were required to assemble and provide 
rigidity to the configuration. In principle there should 
be no problem with the assembly of N-inlet segments, but 
in practice small variations in rod tension and combined 
machining tolerances causes the system to behave like 
the connected linkages of a snake. The configurations 
were assembled on a flat bench and hoisted into position. 
The multiple surfaces were satisfactorily sealed by 
vacuum-greased mylar gaskets between the flat faces. 
Pressure and flow data were recorded using automatic 
digital equipment (7,8). 

In general, the agreement in mass flow rates as 
metered by the tank venturi and the exhaust orifice is 
±5 percent. The pressures and temperatures are within 
±1 percent, and systematic deviations between pressure 
transducers are less than ±0.5 percent. The working 
fluid is nitrogen, and the reduced inlet stagnation tem- 
perature ranges from 0.7 T r 0 to ambient gas with the 
reduced inlet stagnation pressure to P r Q < 2.5. 

RESULTS AND ANALYSIS 

The theory for flows through sequential apertures 
(e.g., orifices) is not well developed. Previous success 
using a combined thermodynamic and choked-flow analysis 
will be applied here (6,7^,£). The governing equations 
are given in the appendix and a description of the theo- 
retical iterative approach used in solving the 
N-sequential-inlet problem follows. 

The process at the ith inlet (Fig. 5) is assumed to 
expand isentropically through that inlet, followed by 
isobaric recovery in the spacer reservoir with (or with- 
out) heat addition. The solution is complicated by the 
fact that the pressure ratios across each inlet, c.- = 
p r,i /p r,i-l* are unknown, and stability of the solu- 
tion is not guaranteed. Thus, to solve the problem of 
flow through N inlets (Fig. 5), one must first determine 
the inlet properties for the first inlet at the given 
stagnation conditions. Next assume a pressure ratio 
across the first inlet, c i * P r ,i/P r , 0 ; then cal- 
culate the conditions of the expansion and determine the 
stagnation conditions for the next inlet (including heat 
addition or extraction in the reservoir spacers). Repeat 
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for the next N-l inlets. At the Nth inlet, determine 
the expansion conditions and compare to those calculated 
from the choked flow constraint, including the effects 
of two-phase flows; if the flow rates and the pressure 
ratios are within a convergence range, they are said to 
be computed for the prescribed configuration of apertures 
( e.g. , orifices) and any addition or reduction of heat 
in the spacer reservoirs. 

In essence, the constraints form the basis of a 
variational approach in which a solution with the least 
increase in entropy is formulated with real fluid pro- 
perties determined using the GASP code (10) . 

Choked Flow Rate 

In theory, similarity exists via the conservation 
equations (see also ( 1-5) ) . For the N-sequential-inlet 
problem, nonsimilar terms are lumped into the flow co- 
efficients. Losses are also lumped into the flow co- 
efficients, with the assumption that they remain nearly 
constant for all stages. This makes the problem tract- 
able, but masks the complexity. 

As pointed out earlier, similarity relationships 
between mass flux and reduced pressure have been demon- 
strated for N < 4. We must now establish form similarity 
for N > 4; this can be done as follows. As seen from 
Fig. 6 by comparing data for N - 7, 10, 15, and 20, the 
reduced mass flux locii are form similar with reduced 
inlet stagnation pressure P r 0 for selected varia- 
tions in the parameter,' reduced inlet stagnation tempera- 
ture T r>0 . This leads to normalizing the flow rates 
with respect to that of the first inlet as a function of 
N and the inlet stagnation conditions upstream of the 
first inlet (T r#0 , P r , 0 ). 


The function f is nearly exponential with N and 
is weakly dependent on inlet conditions (P r 0 ,T r 0 ). For 
many cases, the exponent may be considered a constant 
ranging from 0.35 for N < 3 to 0.45 for N > 15. Thus 
one may write: 

f(N) - N- 4 (3) 

A slight improvement comes from a more complex quadratic 
relation. However, one can use the data figures 
directly. 

Superimposed on the data of Fig. 7 are the theo- 
retical locii calculated using the iterative procedures 
of the appendix and outlined previously. As 
can be seen, the data and the analysis are, in most 
cases, in reasonably good agreement. Moreover, the 
change of slope corresponding to the exponent variation 
(Eq. (3)) is also predicted by theory. As the results 
of the analysis have previously been shown to possess 
form similarity (valid for even such diverse apertures 
as the orifice and Borda types), it appears that extend- 
ing the similarity hypothesis to the case of N sequen- 
tial inlets is justified. 

With Eqs. (1) to (3), one can predict the mass 
flux for N sequential inlets over a large range in inlet 
temperature and pressure and possibly for a variety of 
fluids when either G r j or Cf is known. 


Pressure Profiles 


' fi r,, ■ 3 ». %,0 P r.o> 


G r,» ' G r,l ■ ' (»• P r,.> 


( 1 ) 
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where 

G r,l ’ C f G r,v 


The function g represents a loss coefficient for the 
configuration, and if the combined flow coefficient and 
friction factor losses for each inlet reservoir are 
similar, f applies. Even though much theoretical work 
has been done on flows through apertures, G r i (the 
flow through the first inlet) usually result^ from a 
calibration facility where it is related in terms of 
flow through the classic venturi G r v using a flow 
coefficient Cf. For the venturi, G’ v can be pre- 
dicted using the two-phase choked-fidw theory (9) and 
with the similarity theory using corresponding states 
( 1-5) . Consequently if the function f were known and 
G r i was either estimated from theory or obtained by 
calibration or through other sources, then one could 
predict the mass flux through N sequential inlets over 
a wide range of temperatures and pressures for a variety 
of fluids. 

These observations lead to the formulation of 
Fig. 7, which illustrates the similarity between reduced 
mass flux with the number of inlets for reduced inlet 
stagnation pressures of 0.7, 1.0, and 1.5 at selected 
values of reduced inlet stagnation temperature. 

Although the levels change with pressure, the locii are 
form similar at these pressures. 


Typical axial pressure profiles for the 
N-sequenti a 1-orifice-in let configurations with inlet 
stagnation conditions (T r 0 ,P r 0 ) are given in Figs. 8 
and 9. In both figures, f>art (a) represents gas flow 
and part (b) represents liquid flows. The abscissa 
represents the number of orifices in the configuration 
or the axial distance (e.g., 155 cm for N » 10). In 
Fig. 8 the spacer-reservoir pressures are connected by 
dashed lines and the orifice pressures are connected by 
solid lines. The pressure profiles between the orifice 
and the spacer reservoir (Fig. 9) can be inferred from 
the more detailed experimental work of Refs. 6 and 7, 
and although such detail is not part of this experiment, 
the profiles can be represented by a series of 'horse- 
shoes' (Fig. 9). It should be immediately apparent 
that the locii of Figs. 8 and 9 for the gas and liquid 
cases differ significantly. In Figs. 8(b) and 9(b), 
the orifice locii vary almost linearly (slightly 
concave) with axial position for each value of N 
tested, whereas in Figs. 8(a) and 9(a), they are convex 
or parabolic and are typical friction loss character- 
istics. The spacer-reservoir locii appear parabolic in 
either case, but are essentially linear for liquid 
flows. Such pressure profiles are not only character- 
istic of N sequential apertures, but also of flows 
through high L/D tubes (12) and shaft seals of high- 
performance turbomachines (13,14). 

The axial pressure profiles for inlet stagnation 
conditions (T r 0 , P r p) appear to be similar for large 
N when normalized in terms of the number of sequential 
inlets and the inlet stagnation pressure. 

(Pr^o * Pp(x=0) ( 4 ) 



X 

X choke 


N j 

nr 

s 


N i 


(5) 


where P r - P/P c , X the axial distance, X <-hoke 
the distance to the point of choked flow (assumed to be 
in the exit orifice), and Nj represents a segment of 
an N-sequential-orifice-inlet configuration (Fig. 3). 
Applying Eqs. (4) and (5) to the data of Fig. 8 
produces the two distinct locii for all values of N 
shown in Fig. 10. At this time, the limitations of Eq. 
(5) are not clear; e.g., the limits of N and that all 
the distances between Nj and Nj_i need be the same. 

Since the liquid and gas profiles represent 


T r 0 * 0.7 and T r 0 * 2.3, respectively, it is quite 
apparent that the’shapes of the locii of Fig. 10 depend 
on T r q. From earlier work, ( 13-15 ) it is known that 
these ’loci i are nearly linear for Y r Q < 1; parabolic 
for T r0 > 1; and experience a transition from linear 
to the parabolic near T-. 0 * 1. However, at 16,200 
L/D (13), the pressure I6cii are a little more concave 
than either here or for the shaft seals (13-14). Such 
evidence indicates that for very large N, the liquid 
profiles will also be more concave. Such corroborating 
evidence also extends the hypothesis of the similarity 
relationship between N segmented rough tubes and N 
sequential inlets. 

It is well known that mass flux in choked flows is 
more readily predictable than are the pressure profiles. 
Yet in terms of the normalized coordinates of Eqs. (4) 
and (5), the pressures do appear to be predictable 
directly from the data or within the limitations of the 
theory as 


P r,j / 


r,o 



( 6 ) 
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where q is the calculated pressure ratio across the 
ith aperture. The theoretical locii (Fig. 10) were 
calculated for N * 10. Thus form similarity for N 
sequential inlets appears to exist for both choked-mass 
flow rates and associated pressure profiles. 

SUMMARY 


We have presented critical mass flux and pressure 
profile data for N-sequential-orifice-inlet configura- 
tions where N * 20, 15, 10, or 7. The orifices were 
uniformly spaced at 15.5 cm, and the working fluid was 
nitrogen. 

The data correlated well over a large range of 
reduced inlet temperatures from 0.7 T r ^ to ambient 
and reduced pressures from near saturation conditions 
to P r#0 =* 2. Data are also in good agreement with pre- 
vious’studies of one to four inlets. 


The agreement between experiment and analysis is 
good for the liquid and gas critical mass flux cases, 
but inconclusive for the near-thermodynamic critical 
regions. 

Pressure profiles normalized in terms of inlet 
stagnation pressure and axial distance to choke appear 
independent of the number of inlets. The pressure pro- 
file shape changes from slightly concave to convex 
(parabolic friction loss form) with increased inlet 
stagnation temperature. 

It does appear that a form similarity hypothesis 
can be extended to the case of N sequential inlets and 
used by the designer to determine both the choked-mass 
flow rates and a description of the associated pressure 
profiles. 


Pressure profiles normalized in terms of inlet 
stagnation pressure and axial distance to choke appear 
independent of the number of inlets. The pressure pro- 
file shape changes from slightly concave to convex 
(parabolic friction loss form) with increased inlet 
stagnation temperature. 

It does appear that a form similarity hypothesis 
can be extended to the case of N sequential inlets and 
used by the designer to determine both the choked-mass 
flow rates and a description of the associated pressure 
profiles. 

SYMBOLS 


Cf flow coefficient 

D diameter of orifice or tube 

f function defined in Eq. (1) 

G mass flux 


G* 

9 
L 

t 

N 
P 
T 
X 

z 

l/v 

c 

Subscripts: 

c 

i,j,M,N 

0 
r 
v 

1 



o 

flow normalizing parameter, 6010 g/cm -s for 
nitrogen 

function defined in Eq. (1) 

tube length 

orifice length 

number of inlets 

pressure 

temperature 

axial location 

compressibility 
density 
pressure ratio 

thermodynamic critical 

ith, jth, Mth, Nth sequential inlet 

stagnation 

reduced by normalizing parameter 
venturi 

case for N » 1, the single inlet, or unit 
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Appendix 


Example of Analysis— An Analysis of Sequential Inlet Flow* 


The flow process is neither steady, adiabatic, nor 
frictionless, as described in the text. However, if we 

u i 

Pi 

Pi 

u i+l 

Pi+l 

r i p i+i 


simply ignore these problems, the governing 
equations for the 1 th inlet, see sketch and figure 5, 
can be written as (ref. 5): 

v v m L — L ! — L j L. 





i : 


Continuity ^IsO 

dx t 

(Al) 

Jetting effects in sequential Borda inlets 


Momentum dpu/ttj _ dp_ ^ 

dxj dxi 

Energy JL P u,H 0 *0 

(A2) 

(A3) 

/oa.O-(^Xa.O 

J o 



where 



► 

(AS) 

H 0 *//+ jU[Ui 

(A4) 

haa '° m 



State p=p(ft,H,x a ) 

(A5) 

Isentropic restraint 



Similarity (extended van der Waals 
one-fluid model for corresponding states) 


TdSi =dH— =o 
P 

i*l, n 

(A9) 

Z m (T,V,x a )-Z 0 (%,p 


Isobaric restraint 



Im(KT.x a )=Sgo(?,p 


Si+ 1 *S(Pi,Ho) 

i- 1,« 

(A10) 

+ K £ x a In x a - In h ) 

(A6) 

Choking constraint 



where 


oL ,-Af) - - 1 

max 'dp' e y2 I 

\ P ° Vdp 
] P. 

t 

(All) 






■A 2 HgXaXtff a /3, O^ 1 a/3,0 


It appears that at the 30-diameter separation 
choking can occur at the entrance or exit of the last 
(/an) of the sequential inlets. At the 0.8-diameter 
separation choking can occur at the inlet of the first 
(/ * 1 ) inlet of the sequential inlets or at the exit of the 
last (/an) inlet. The details of the choking con- 
straint, at a fixed position, for flows with change of 
phase are quite complex and, although used herein, 
will not be repeated at this time; see reference 4 for 
further details. 


♦Reprinted from NASA TP-1792, May 1981, with changes in 
in ref. and fig. numbers. 
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Furthermore, if one isolates an inlet (one section of 
fig. 3 or 5) and forms it into a blackbox, the 
conservation equations yield the isentropic expansion 
restraint 


G 2 


td±i [ Pl d ± 


(A16) 


dHo-dH+d- 


dp+pd-^ *0 


or 


Again, consider two cases: (a) when Aj> >A i+l , 
case (1) above results; and (b) when 
M 1 + 1 P /+1 -P/A,), it leads to case (2). 

(A12) (4) In jetting data (refs. 3 4, and 13) p i+i >p h 

and it follows that u ,+ \ <ui and that the system 
functions as a diffuser. 

The same results are achieved with enthalpy as the 
independent variable; from the energy equation it 
follows that 


dp 

P 


—dH=~0 


«? +1 -^-2(H 0 -// <+ j)-2f' ^ (A17) 

J P/+ 1 P 


Thus wc need only consider an isentropic process, 
and from the momentum equation we find 




(A13) 


Consider now the following cases: 

(1) When Mj+i > >Uj or u,_o, the mass flux 
across the i 11 * inlet becomes 


G 2 *(P/+i«i+i) 2 *2p? +1 P ^ 

J Pi*i p 


(A14) 


(2) When u i+l -u, or u,-u /+1 , p i+l -pi and 
jetting can occur, provided the jet has been 
established elsewhere in the system and 


The enthalpy relation is used for the preliminary 
calculations herein. 

In this report we use the simplified form of van der 
Waals corresponding-states principle since we only 
have data for fluid nitrogen. 


^aa,0 * 1 


*\ 


<t>aa,0 ■ 1 
f=T R = 


L 

T c 



> 


h = v *=V c J 


(A18) 


and 





•A' 3 * C„= £ 

R G* 


(3) For plug flow continuity, w=p,u;.4, = 
P/+i“/+iA,+i, the mass flux becomes 


All thermophysical properties were calculated by 
using the computer code GASP (ref. 10). 
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Figure 1. - Schematic of flow facility modified for operations with N sequential inlets. 






Figure 5. - Schematic H-S diagram of the N-sequential-inlet flow 
process. 
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Figure 6. - Reduced mass flux as a function of reduced inlet stagnation pressure at selected reduced inlet stagnation temperatures, 
T r 0 for N sequential orifice inlets with 15.2-cm spacer reservoirs. 
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PRESSURE, MPa 


NUMBER OF 
ORIFICES 




(b) Liquid flow. 


Figure 8. - Average connected-orifice and spacer- reservoir 
axial pressure profiles for 7, 10, 15, and 20 sequential 
orifice inlets for gas and liquid flows. 
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SPACER-RESERVOIR 

AVERAGE 




(a) Gas flow. 


Kg 

Ifc 




'SN. S 






\V 


5 10 15 2 

AXIAL DISTANCE, EQUIVALENT TO (N-l) x 15. 5 cm 

(b) Liquid flow. 

Figure 9. - A comparison of axial pressure profiles based on detailed measurements of Ref. 7 and average connected-orifice and 
spacer- reservoir profiles for 7, 10, 15 and 20 sequential orifice inlets for gas and liquid flows. 





NORMALIZED AXIAL POSITION, (X/X choke = Nj/TNj) 

Figure 10. - Normalized axial pressure profiles for N sequential 
orifice inlets for gas and liquid flows. 
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